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Abstract 

We calculate Higher Order corrections to the resolved component of the 
electroproduction cross section of large- E"^ hadrons. The parton distributions 
in the virtual photon are studied in detail and a NLO parametrization of 
the latter is proposed. The contribution of the resolved component to the 
forward production of large- i?^ hadrons is calculated and its connection with 
the BFKL cross section is discussed. 



1 Introduction 



The electroproduction cross section of \arge-E± hadrons can be split up in two 
parts. One of them describes the reaction in which the initial virtual photon takes 
part directly in the hard scattering process ; it is called the direct part. But the 
photon can also act as a composite object which is a source of collinear partons 
which will take part in the hard subprocess ; this mechanism is usually refered to 
as the resolved process and defines the parton distributions in the virtual photon 
which have the feature of being proportional to £n Ej_/Q'^ in the asympototic region 
where Ej_ ^ (virtually is the absolute value of the photon). 

This distinction between direct and resolved component parts is especially useful 
in photoproduction reactions in which a quasi-real photon is present in the initial 
state (for a review, see ref. PP). In this case the parton distributions in the real 
photon are proportional to in E\j h?QQj^ and can be quite large. The interest in 
these real distributions dates from the pioneering work by Witten [2j who showed 
that their asymptotic behavior can be completely calculated in perturbative QCD, 
a result which opened the way to interesting tests of the theory. Nevertheless, when 
-^i/^QCD decreases, the importance of the non perturbative contributions grows 
and we return to a situation similar to that of the proton structure functions for 
which non perturbative inputs are necessary. 

The situation is clearer when the initial photon is not real, but has a virtuality 
much larger than Aq^^. In this case the non perturbative contributions (for instance 
that of the Vector Meson Dominance type) are suppressed by powers of and we are 
back in the realm of perturbative QCD. The magnitude of the virtual distributions 
is smaller than that of the real distributions. Nonetheless, they are observable 
and dedicated experiments have studied the virtual parton distributions in e^e~ 
collisions |H1 E] and in the electroproduction of large E^ jets El E! and hadrons 
[HI E] ■ These studies acquire a quantitative status when data are compared with 
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theoretical predictions calculated beyond the Leading Logarithm approximation |10|, 
El El El E]- It is the aim of this paper to establish such NLO expressions for the 
resolved component of the electroproduction of large- ii^^ hadrons. We studied the 
corresponding direct component in ref. [T^ . 

This work puts the theoretical predictions on a firmer ground since the full cross 
section formed by the direct and the resolved component parts is now calculated at 
the NLO approximation. In ref. ^Sl we founded predictions for the leptoproduction 
of forward large-i?_L hadrons on a NLO calculation of the direct term only. Then 
we observed that the resolved component, calculated at the lowest order, was not 
negligible. Here we pursue this study of the forward production now including the 
HO corrections to the resolved part. This allows us to refine our predictions and 
our comparisons with the BFKL-type cross section which should constitute a non 
negligible part of the forward cross section ^t^l IT7j. 

In the next section we gather kinematical definitions and general expressions 
concerning the resolved cross section, including a discussion of the kinematical do- 
main in which such a resolved component can be defined. Section 3 is devoted to the 
general structure of the NLO corrections and the issue of the factorization scheme. 
In Section 4 we propose a parametrization of the NLO parton distributions in the 
virtual photon, finally, we consider some numerical applications in Section 5. 

2 The resolved component 

In this section we present the kinematical definitions and the general expressions 
necessary for the study of the resolved component. This determines the frame in 
which the HO calculation described in the next section, will be performed. The 
cross section of the reaction e(£) + p{P) e{i') + h{P4) + X, 




(1) 
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is written in terms of the leptonic tensor i'"' = 2{i>'i"' + i^i'^" - g'"'{i ■£' -ml)) and 
of the hadronic tensor T^^, which describes the photon-proton coUision. We define 
the photon variables = —q^ = —{i — i')'^ and y = ^^ff^ = pif = Q'^/i^BjS) in a 
frame in which has no transverse component (we neglect the proton mass and 
is positive (HERA convention)). S is given hj S = {P + i)"^ and Xg^ has the 
usual definition = Q^/2P ■ q ; (f is the photon azimuthal angle. The differential 
phase space of the final hadrons is given by (a sum over the number of final hadrons 
is understood in ((H)) 

dPS = i2n)H' {q + P-tpi^ n • (2) 

The hadronic tensor can be calculated as a convolution between the partonic ten- 
sor tfj_i, which describes the interaction between the virtual photon and the parton of 
the proton, and the parton distribution in the proton Ga{x, M). The fragmentation 
of the final parton which produces a large- hadron is described by the fragmen- 
tation function D^{z, Mp). These distributions depend on the factorization scales 
M and Mp, 

I T^^dPS = Y.j —Ga{x, M) I dz Dl{z, Mp)tf^ ■ dps (3) 

a,b 

where dps is the phase space element of the partons produced in the hard photon- 
parton collision. From expressions (j21) and Q, we obtain 



da E±4 « /■ , ^ / Tt^\ f dz 



J dxG,{x,M) J -^D'^{z,M, 



d(pdQ'^dydE±4dri4 271 2it 

27r (47r)2 2xS q^ ^ ^ 

where the phase space dps' no longer contains parton 4 which fragments into h{Pi). 
{rj4 is the pseudo-rapidity of the observed hadron). 
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It is useful to give a more explicit form to the tensor product in the 7* — p frame 
by defining the transverse polarization vectors Si = (0, 1, 0, 0), 62 = (0, 0, 1, 0) and 
the scalar polarization vector = — i=(g^, 0, 0, with = 0, 0, q^) the virtual 
photon momentum 



+^^i.jQ' tsi + Q'^^^-p^tss , (5) 
y ^ 

the transverse momentum of the initial lepton being along the x-axis. 

In the limit and after azimuthal averaging over ^94 we recover the unin- 

tegrated Weizsacker- Williams expression 



(6) 



2 ^ yQ2 Q4 

with cr_L = ^ (tn +^22). 

Actually the limit (jUj) is correct only if lim tgs = 0{Q'^). This is not true if 
an initial collinearity is present in the partonic tensor (light partons are massless) 
which leads to the behavior lim tgs = This point is discussed at the end of 

this section. 

The partonic tensor is given by a perturbative expression in a^. The Born 
contribution is of order 0{as) and corresponds to the QCD Compton subprocess 
7* + g g + q and the fusion process 7* + 5' — > q + q- Higher Order 0{a1) corrections 
to the Born cross section have been calculated in ref. ^Hl- In the course of these 
HO calculations a resolved component appears, corresponding to subprocesses in 
which the virtual photon creates a coUinear q-q pair ; the quark or the antiquark 
subsequently interacts with a parton of the proton. 

Let us study this contribution in detail by considering the simple model illus- 
trated by the gauge invariant set of Feynman graphs displayed in Fig. 1. The 
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neutral parton of momentum p is off-shell and is part of a hard process also involv- 
ing a parton of the proton. The final parton of momentum p4 fragments into the 
observed large- hadron of transverse energy E±4. All the results described below 
can easily be obtained from the expressions given in appendix 1. 




Figure 1: Feynman graphs leading to a resolved contribution. 

The cross section corresponding to the graphs of Fig. 1 has double and single 
poles in = {q — k'Y. The interference term between graphs (a) and (b) has a single 
pole which leads to an expression proportional to ^n^^, after integration over k^^. 
However, a prefactor (f is present in all tensor components Iab (a, b= S", 1,2). As 
a result, these components have no singularities when tends to zero. This well- 
known behaviour is due to current conservation (for the components involving a 
scalar photon) and to the fact that interference terms are not singular for transverse 
photons. Therefore, let us concentrate on the square of graph (a) and start with the 
transverse component which has the expression (after integration over the azimuthal 
angle ^pk') 

t.. = -3e}/d.[(l-.r + .=]|^^^ -{..--eyj!^ (, = 1,2) (7) 

where //(O) is the hard subprocess amplitude, here representing the process A;-|-p — > 
P4, in which we have set k"^ and k\ equal to zero. The upper limit of the k"^- 
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integration indicates the scale at which the coUinear approximation used in (|7j) by 
setting k"^ and fc^ equal to zero is no longer valid, using formulae (jH) and and 
after integration over k"^. The contraction with the leptonic tensor leads to 



y Q\ 

In-^ -tnz- £n{l - z) - l\ a , 



difdQ'^dy 2tt 

pi 
Q 



Jdz^ 3e}[{l - zr + z^] 



where we define a = ^^^"^ (we have not written the contribution of order 0{Q'^/p]_^)). 
This expression is the lowest order resolved cross section and is exactly the one which 
is obtained in the course of the calculation of HO corrections to the direct Born terms 
|15j . The Weiszacker- Williams distributions of the virtual photon in the initial elec- 
tron and the quark distribution in the virtual photon are universal as they do not 
depend on the particular hard process described by cross-section a. Expression (jHl) 
is the starting point of this paper. Indeed, when in^^ is large, one cannot con- 
tent oneself with this approximation and corrections of the type (^in^^ with 
k = n,n + l must be calculated and resummed. These corrections modify expression 
(jH| at the Leading Order {k = n) and at the Next-to-Leading Order [k = n + 1) 
approximation. 

In order to avoid double counting, expression (jS)) must be subtracted from the 
NLO direct cross section. Actually the exact expression to be subtracted is a matter 
of factorization scheme. We define the resolved component by 



27rda''^' a 



difdQ'^dy 2tt 



l + {l-yY 1 2miy 



r2 



dz — 3ei 

271 ^ 



zf + z' 



Q 



y 

where we introduce the factorization scale with = 0{E±4). After subtrac- 
tion, the part of (jH)) left in the direct HO corrections is obtained from (jH)) by the 

2 2 

substitution in^^ — >■ in^j^. We call this factorization scheme the virtual factorisa- 
tion scheme. This is a natural scheme in virtual photoproduction in which all the 
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In Q^-terms are resummed in the parton distributions. Then the total NLO cross 
section is given by the sum of the subtracted direct cross section and of the resolved 
cross section calculated at NLO at the scale M^. The variations of the resolved cross 
section with are partly compensated by the In terms, these remain in the 
direct cross section so that the total NLO cross section exhibits a smaller sensitivity 
to than the LO cross section. 

Of course this procedure is useful as long as ^ Q^. Actually the coUinear 
approximation used in is valid if Q"^ k\ < p\^, which allows us to put k"^ = 
in the hard cross section. When p'^^ < Q^, this upper limit is incorrect. Let us 
rewrite the fc^-integral in ^ in terms of fc^ 

/ ^^(^i) = / ,2 fn V ^(^1) • (10) 

J-q2z\k^\ Jo kj_ - q^z{l - z) 

This integral is sensitive to the dependence on kj_ of the 2 — » 2 subprocess cross 

section 5(^5.) which behaves approximately like 0{j^^-^^^-^). This behavior shows 

that no collinear logarithmic terms (coming from the denominator k\ + Q'^z{l — z)) 

are present when p5_4 ^ Q^- Therefore, the resolved component must be proportional 

to the result of /c5.-integration in which the upper hmit p5_4 is replaced by +p\a- 

For p5_4 ^ we have the case already discussed and for p\^ <^ Q^, there is no 

resolved component. 

As a consequence it is more appropriate to define the factorization scale 

m2 = g2 + c'^^El, (11) 

(i?_L4 is the transverse energy of the observed hadron) which has the following correct 
properties. 1) It does not depend on kinematical variables internal to the subprocess 
which may lead to incorrect results when HO corrections are calculated [19, ; 2) The 
resolved component calculated at vanishes when ^ ; 3) again we find 
the conventional factorization scale ~ C^E±4 when Ej_^ ^ Q^, Cy being an 
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arbitrary constant of order 1. 

Let us finish this section by discussing the tensor components tis and tgs which 
come from the square of graph (a) in Fig. 1. The components t^g behave hke 

2 

y/q^in^^ and have no singularity at the hmit —* 0. On the contrary tgg has a 
constant behavior when —>■ 

or 

a result which leads to the scalar cross section 

— — , = — dz — SCf 4z 1 — z)a . 13 

In going from (fT2|) to (fT3|) . we dropped the —q^z{l — z)/p\^^ term which depends, 
through px4, on the detailed kinematics of the subprocess. 

We observe that t^s has a "constant" behavior when — due to the double 
pole of the cross section. Actually the limit — ^ corresponds to a non perturbative 
region for the /c^-integration. If instead of Ifc^ml = —q^z we set |A;^j„| ~ ^qcd^ 'we 
would obtain a vanishing cross section when g^ — > 0. A similar result is obtained 
if we consider massive quarks (with — = —q^z + fz^)- Therefore, for a 

physical process and a real photon, there is no tss contribution, as can be expected. 

However, let us notice that for small values of ~ ^qcd^ ^^e resolved cross 
sections, as defined in (jH)) and (fT^ . strongly depend on the way the fc^-integral is reg- 
ularized, different lower bounds produce different 2;-dependence, and thus different 
physical results even when Ej_^/Q'^ is large. This paradox is however solved by the 
HO correction to parton distributions in the photon discussed in the next section. 
There we shall see that the NLO parton distributions contain a term that cancels 
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the unwanted 2;-dependent contribution, up to a vanishing term when Ej_^/Q'^ tends 
to infinity. Actually this result is true for all z-dependent terms of coUinear ori- 
gin (related to the lower limit of the ^^-integration) present in (jHJ and As a 
consequence the scalar cross section (fT^ will be cancelled. 

3 NLO corrections 

In section 2 we defined the resolved component of the transverse cross section 
{i = 1,2). (Here is defined as the Born amplitude squared divided by the fiux 
factor z). 

tu = ^e}ldzPi^\z)in^a- (14) 

where the factorization scale is given by (jllj) and Pg^\z) = 3[2;^+(l— z)^] (for 1 quark 
species). Expression (|TH) contains the lowest order ((9(q:°)) parton distributions in 
the virtual photon 

<l{z,M^^,Q^) = ^e}Pi^Hz)in^. (15) 

The Born cross section describes the scattering between a quark of the virtual 
photon and a parton of the proton producing two large-p_|_ partons in the final state. 

Leading Logarithm (LL) corrections, corresponding to the emission of coUinear 
gluons by the initial quark, can be obtained by solving the following inhomogeneous 
DGLAP equation fH^ ITH] (we only reproduce the evolution equation for the Non 

Nf 

Singlet (NS) quark distribution q^^^ = Qf + Qf ~ (<?/ + Qf)/^/) with < ej- >= 

(16) 
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where Pg(°)(z) = Cf{{1 + z^)/{1-z))+. The lowest order expression (fT^ is solution of 
such an equation when P^^ is set equal to zero. The solution of (fTH|l for the moments 
gf ^(M2, Q^,n) = dz q^^{M^, z) is given by {d{n) = 2Pi^\n)/f3o and f3o 



is the lowest order coefficient of the /5-function expansion g^^?°2s = P{as) = —^Po 



n 



(17) 

with the boundary condition q^^'^^{Q^,Q^,n) = 0. The Leading Logarithm NS 
expression for the resolved cross section is now given by 

e = g^""®?^ (18) 

which is expression ()14|) in which the lowest order parton distribution is replaced by 
the LL solution (jl7p . 

The next step is to look for a Next to Leading Order (NLO) expression for tu, 
which requires the calculations of HO corrections to both q^^ and . Indeed the 
structure of these HO corrections is the following. The hard cross section has the 
expression {a^ is of order 0{al)) 

whereas the parton distributions behave, in the asymptotic domain M^/Q"^ ^ 1, 
hke 

It is clear from (fT^ and (j^n|) that a NLO expression for tu can only be obtained by 
calculating both a'j.B and b. 
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3.1 The hard resolved cross section at NLO 

The calculation of the HO corrections to the hard resolved cross section is the 
simpler part of the NLO program, since these HO are the same in real and virtual 
photoproduction reactions, providing we work in the same factorization scheme. 
Therefore, we can borrow the results of ref. [20! obtained for the real photoproduc- 
tion of large E_i_ hadrons. 

Let us elaborate this point by first studying the resolved Born term. In the real 
case, instead of ((Zj) we obtain the following expression 

t^^- ^^6efJ dz[Z +[l Z) r(l- 25)70 {klV^^^^ ^ ^ 

in which we use the dimensional regularization and n = 4 — 26. The expression 
between the square brackets is the n-dimensional DGLAP branching function ; the 
factor (47r)^r(l — e)/{r{l — 26){1 — e)) comes from the azimuthal integration and 
the n-dimensional photon spin average. After integration over kj_, we obtain {1/e = 
I + in An- ^e) 

(fe^lnV^) - P°r,{z) + 3) S'' (22) 

where the limit, when e tends to zero, of the n-dimension Born cross section af is 
simply of expression (fT^ . This expression is identical to that obtained in the 
calculation of the HO corrections to the real direct term. 

At this point, if we subtract the term proportional to ^— ^ + ^n,—^ from (j22j) . 
which defines the MS factorization scheme, we obtain a direct HO subtracted con- 
tribution different from the one found in the virtual case (cf. expressions (jHI) and 
0). However, as we shall see in the next subsection, this scheme dependence is 
compensated by the NLO corrections to the parton distributions. 
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Now let us go one step further and consider 0{as) corrections to the resolved 
expression 0. These HO corrections are the same in the real and in the virtual 
case, with the exception of collinear contributions coming from the branching 7* — » 
q + q + g and containing [in p'x^/Q^)"' {n = 1, 2) terms. These logarithmic terms can 
be factorized and resummed at the NLO approximation with the result (we consider 
only the Non Singlet case) 



where ® indicates convolutions in the longitudinal variable. The factor 1 in the 
parenthesis corrsponds to the Born contribution Q. The term ^hgg (g) is the 
collinear HO correction calculated in the virtual factorization scheme (resummation 
of all the in terms in the parton distribution function with the boundary con- 
dition q^^^{Q'^,Q'^) = 0). However the physical (direct + resolved) cross section 
is factorization scheme invariant and can be written in terms of the MS quantities 
q!^^'^ and hqq{z). As a result we can use the HO correction calculated in ref. [201 in 
the MS scheme if we also use parton distributions (and a direct term) calculated in 
the same scheme. 

The authors of ref. ^21 ^1 ^^^'^ worked in the MS scheme in their study 
of the elect ropro duct ion of large-p^ jets, and they established the expression which 
must be subtracted from the virtual direct term in order to obtain the MS direct 
term. We comment on their results at the end of section 3.2. 

3.2 The virtual parton distributions at NLO 

In order to delimit the problem of the Factorization Scheme (FS) in the virtual 
parton distributions, we study the simple case of the DIS on a virtual photon and 
we consider the n-moment of the structure function J-"] = F2{x, K^, Q'^)/x in which 

= Ig^l is the virtuality of the target photon, = \k'^\ is the virtuality of 
the probe photon and x the Bjorken variable. To make the connection with the 




(23) 
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transverse cross section defined in (|7j), JF^ is defined by an average over the transverse 
spin of the target photon only. To simplify the discussion we only consider the Non 
Singlet contribution. T2 is the sum of a resolved part and a direct part (we drop 
the indices n) 

n{K\ Q') = C2,,MK')) . q!^'{K', Q) + C^^^{as{K')) . (24) 
In ()24|) is proportional to ^fl^f^ < >] ; we drop this factor which is 

/=! 

useless in the present discussion. The direct part, C^^, and the resolved hard cross 
section C2,q (the Wilson coefficient) are expansions in agi^K"^). All the Cn depen- 
dent terms are collected in the virtual quark distribution g^"^ with the boundary 
condition q^^ {Q"^ ^Q"^) = 0. This defines the virtual factorization scheme already 
mentioned in section 2. In fact ()24|1 is the final result obtained by Uematsu and 
Walsh ^U] in their study of the virtual photon structure function, using the OPE 
and the MS factorization scheme as a starting point. The distribution g^"^ verifies 
the inhomogeneous DGLAP equation (fTT)|l in which the lowest order branching func- 
tion Pj°) and Pg^g^ must be replaced by the all order functions P^'^(as(M^)) and 
P^'^(as(M^)) = "''(^-^^ ^ Pqq + ■ ■ ■ which are expansions in as(M^) and depend on 
the factorization scheme. The solution of can be written (from now on we drop 
the index NS) 



a rMK^) da' P,,{a') |^p(„. 



Til being a physical observable, must be FS scheme invariant and cancellation 
must exist in ()24|) between the various scheme dependent contributions. Let us first 
note that C^^-^iK'^^ is FS scheme invariant because 

T''^{Q\Q^) = C<,^,{aXQ''y) . (26) 
To study the scheme dependence of g^, let us start from expression (j25p and 
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define a new DGLAP branching function P by 



P^P-5P 



where SP is an arbitrary expansion in as starting at order 0{al) 



27r 



with given by 



(27) 



(28) 



J a 



^s{K^) da' 



J a' /3(a") 



(29) 



,(Q2) l3{a') 

We see that the variation 5P can be absorbed in C2,q (the hard resolved subprocess) 
and Pg-y, thus defining new expansions in ag, C2,q{oisiK'^)) and Pq^{a'), whereas 
is kept unchanged 



Let us now study the effects of modifying Pq^ 



(30) 



P — P — SP 

J- q-y — ?7 "-'-'97 



with the arbitrary series 5Pq^ starting at order 0{as) 



a ns{K^) da 



5Pn^ eJc' f> 



97 



where the parton distribution q is calculated in the bar-scheme 



(31) 



(32) 



/.^9 « rM^^) da'— r- 



) dc"- 



(33) 



Finally for we obtain the expression 
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where 



C2,,{asiK')) = C2,y{asiK'))-C2,,MK'))— / —6P,, eL' (35) 

/TT JO p 



and 



<(i^^Q^) = ^/ ^5P,,eL' . (36) 

Therefore in the new factorization scheme (the bar-scheme), the structure of the 
expression for is the same as in the original scheme, but the parton distribution 
does not vanish at = since q^iQ'^iQ'^) is different from zero. Therefore, by 
going from the virtual FS to the bar-scheme, we find the boundary condition that 
the bar- distribution must verify. By rewriting ()36p as 

q^{K\ Q2) = eJ-w^) - / -^^n. eJ"' , (37) 

we see that q^^{K'^,Q'^) verifies the homogeneous DGLAP equation and that the 
boundary condition is given, at the lowest order (P = ^Pqq\ /?(««) = —^Po and 

«,W.Q^) = -^^. (38) 

The bar-scheme can be any scheme, but it is convenient to work in the MS 
factorization scheme in which the two-loop branching functions and P'^^ are 
known. Moreover, in the electroproduction of large-p_L hadrons we also know the 
NLO resolved subprocess cross section (the equivalent of 6*2,5) calculated in the MS 
scheme in ref. j2^. It is easy to obtain SPjj}^^ from expression ()35|) written at the 
lowest order on as 
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C2,7 - C2r, = — ■ (39) 

6*2,7 direct term [22] and 6*2,7 ^'^^ virtual-scheme direct^ term ^U] 

6*2,7 = |8x(l - x) - 2 + (x^ + (1 - x)2) , (40) 

which leads to 

q^{Q\Q') = -^6{{x' + {l-xr)in[x{l-x)] + l} . (41) 

Let us finish this section by going back to the direct cross section of large-pj_ 
hadron electroproduction. The MS boundary condition can be obtained by compar- 
ing expression (jH)) calculated in the virtual case and expression corresponding 
to the real case. We see that subtracting (z^ + (1 — ^Y)^-''^ z{i-z) ~ ^ from the virtual 
expression, we find the MS expression 1 — (z^ + (l — z)). The term that we subtracted 
is equal to (PT|) as can be expected^. 

3.3 The scalar parton distribution at HO 



In section 2 we found a scalar resolved contribution ()T3|l to the electroproduction 
cross section corresponding to the scalar distribution Qq^z) = ^3eg[4z(l — z)]. HO 
corrections to this distribution correspond to the Feynman graph 
of Fig. 2 (with an extra gluon in comparison to Fig. 1). Working in the LL approx- 
imation and considering only terms proportional to in K'^/Q'^), we have 



gf (x, K^, Q') = Jdz q^iz) J dz' ^ Pi^\z') 5{zz' - x)in K^/Q^ , (42) 

"'^This direct term corresponds to transversely polarized photons whereas the expression of ref. 
[TU| also contains the scalar contribution. 

^The subtraction term established in ref. ^21 is identical to the one found here except for a 
term proportional to [z^ + (1 — z)'^]£n z. Therefore it does not totally ensure the transformation 
from the virtual scheme to the MS scheme (for instance from C2,-y to C2,-y in the DIS case). 
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Figure 2: A higher order correction to the Feynman graphs of fig. 1. 
OTqf{n,K\Q^) = q^{n)^Pill\n)£nKyQ\ The full LL expression g^(n) = § ql{n) 

k=l 

is easily resummed 

q'{n,K^Q')= , ^^o^H f e^"' (43) 

JosiQ^) P[0! ) ATI 



This solution is similar to expression (j25jl . but with an inhomogeneous branching 
function starting at order 0{as) ; it can be written 



q'{n,K\Q'') = -q^in 



1 - 



(44) 



Adding this contribution to the lowest order one (fTSjl . we see that the latter is 
cancelled and replaced by a contribution which vanishes asymptotically because of 
the factor (as(-^^)/tts(Q^))~^'^''''^^"- Therefore the scalar contribution ()13|1 which is 
target dependent (it depends on the regularization of the fc^-integral as discussed in 
section 2) is cancelled. This mechanism is actually quite general and also valid for 
the transverse component. 

Therefore, in the electroproduction case a full treatment of the scalar cross sec- 
tion amounts to subtract expression (fT^ from the NLO direct cross-section and to 
add the scalar resolved component 
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^^scalar a 2{1 - y) 1 



— J dzq'^{z,MlQ')a (45) 



d(pdQ'^dy 2n y Q 
where we define^ 

g>,M,^Q^) = go^(n)(^^^^j . (46) 

4 NLO parametrization of the virtual photon struc- 
ture function 

Although the works of Uematsu and Walsh JU] and Rossi [25] date back to the 
eighties, interest in the virtual photon structure function grew much later, thanks to 
the HERA experiments. Since then several papers have been published discussing 
the parton distributions in virtual photons at the LL approximation [211 12Z1 12H] and 
NLO approximation [201 EDI E], and emphasizing the possibility to measure them in 
electroproduction experiments [13 CH E2 IHHl OH IHH IHE] ■ In this section we present 
a study of the parton distributions in virtual photons performed in the MS scheme, 
using the results of section 3 on the inputs at = Q^. We choose large enough 
to neglect non perturbative effects, and therefore we do not present results on the 
limit ^ 0. A NLO study has also been done by authors of refs. [21] and [HO] in 
the DIS^ scheme with emphasis the real limit 0. We may note however that 

the solutions of [30] do not fulfil condition (pUj) of section 3. 

The parton distributions are solutions of the same NLO inhomogeneous differ- 
ential equations as in the real case. The NS equation is given by expressions in 
which ^Pq^^ and ^Pgg^ must be replaced by NLO (two loops) DGLAP kernels ; the 
equations of the singlet sector as well as the expressions of the kernels can be found 
in ref . [21] . The only change with respect to the real case is the starting point of the 
evolution, instead of ~ (.5)^ GeV^ in the real case. The boundary condition 
■^This parton distribution in the scalar virtual photon has been studied in ref. [23 EJ- 
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for the distributions are given by expression ()41|) for the quark distributions and 
we have g^{Q'^,Q^) = for the gluon distributions. Another difference from the 
real case is the existence of a scalar contribution. The latter has been discussed in 
section 3.3. 

In section 3 we studied the massless case <^ Q^. However \^Q^ can be smaller 
than the bottom mass mf, and we have to determine what are the relevant boundary 
conditions for the bottom quark distribution in the virtual photon. To find these, 
let us introduce in expression the kinematic corresponding to the case in which 
the photon interacts with a massive quark. 



3 / _ . „ „ , „ N rh in 



We notice that for = 0, we find the massless corrections already given in (|H)) 
which are associated with the virtual factorization scheme. 
For = 0, we find 

l|p(;,)£n% + 6^(1-^)1 (48) 

which is the MS correction (once the term is subtracted). Therefore, in the 

massive case (Q^ = 0) we directly work in the MS-scheme and we do not have to 
modify the factorization scheme as discussed in section 3. Therefore, there exists a 
transition between case ^ and case ^ that we should study in detail. 

Let us start from expression (jTfjl . When < m^, we factorize Pq^ln-^ which 
is the contribution given by the evolution equation starting at the scale m?. The 
rest is given by (without the 1/ z prefactor) 

2 

^57^^f^ + 6^(1 - ^) + MC<{m\ Q\ z) (49) 
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with 



m 



2 



3-Q2^(l 



+ Q^z{l — z) + Q^z{l — z) 

namely the usual masless MS correction and corrections in jrr?. When > m^, 
we factorize Pqjln-Qj- and we obtain 

P11 ^^T^ + Pi-r^^ ,} N + 6;2(1 - ;2) - 3 + MC>(m^ Q'; ;2) (50) 

-2(1 — 

with 



, n „n , ^ ^ — 3-m 



'^^ m2 + g2^(l -z) m2 + z(l - 2)g2 ■ 
We recognize the massless corrections in the virtual scheme and a massive m^/Q^ 
correction. 

However the same massive corrections MC^{m^Q'^\ z) appear in the calculation 

< 

of the inhomogeneous kernel k^^^ as outlined in section 3.3. When we add the resolved 
and the direct contributions, these massive corrections are cancelled and are replaced 

by MC<(m2,g2.^)(^^(jvf^)2/^^(^2))-2Pi°V/3o or by MC>{m'' ,Q''-n){as{M^) / 
tts(Q^))"^^™'^'^" when Q'^ > m?. In the latter case, we still have to add —[Pq^{z)tnz{l- 
z) + 3] to MC>{m?, Q"^] z) to move to the MS scheme with the result 

MC>(m^ Q'- z) = f ^ - 2Tn2\T^^ ^ (^1) 

m^ + Q^z[l — z) m^ + Q^z[l — z) 

which is equal to MC<{m?, Q^; z) when — Q'^. 

Therefore, we can summarize our treatment of the massive quarks in the following 

way. First we assume that the m^/Ej_^ corrections are properly taken into account 

in the direct term (or that they are negligible when £"14 3> m^). Second, we work in 

the massless MS scheme and we take into account mass corrections through the input 

of the quark distributions. Thus we have the following inputs for < < ml 

(up to charge factors) 
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2^ 



u{x, Q 
s{x, Q 



2\ 



c(x, Q^) ~ Pg^ in 



3g2^(l - z) 



ml + Q'^z{l — z) ml + Q'^z{l — z) 

3Q22(1 - z) 



b{x,ml) ~ Pqy \ FrT^f7T~^ ■ (^2) 

mf, + Q^z{l — z) mf^ + Q^z[l — z) 

Whereas for > m^, we have the input 

b(xO^)r.P in '^'^^^ " (53) 

b{x, Q)r^F,, in ^ ^^^^^ ^ ^^^^^ _ . (5dj 

With these inputs we obtain the distributions shown in Fig. 3. We have chosen 
= 8 GeV^ and = 25 which correspond to average values of and = 
+ of the HI experiment jH] . The distributions calculated in the MS scheme 
increase for x going to one as we can see from Fig. 3. This increase is compensated 
however by the behavior of the direct term which contains terms in in{l — z) that 
become negative at large z. We also remark the effect of the massive input (j52p for 
the charm quark distribution. 

We end this section by comparing our results with experimental data obtained 
by the L3 collaboration P| for the structure function 



Flfs = + Fl (54) 

where the indices _L and s refer to polarization of the target photon of virtuality 
(called in ref. |3]). In terms of these components, the usual structure function 
F2 is written (the tensor indices refer to the target current) F2 = — |5'"^(-^2^)o/3 = 
F2 j_ — 1-^2^8- Until now we calculated only the transverse distributions (Fig. 3) ; 
in order to obtain Fj^-j we have to add to F^ j_ the scalar contribution defined 
in expression ()46p for the quark component, in which a gluon distribution is also 
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Figure 3: The parton distributions in the virtual photon for = 8 GeV^ and 
M2 = 25 GeV^. 

generated by the DGLAP evolution equation. All our calculations are done in the 
MS scheme. 

In Fig. 4 we see that our predictions are in reasonable agreement with data at 
low and medium values of x, but they undershoot them at large values of x. Similar 
results have been obtained by the authors of ref. 

5 Numerical results 

We now turn to a phenomenological study of the Deep Inelastic Production of 
laTge-E± hadrons. We concentrate mainly on the resolved contribution studied in 
this paper and consider the HI data |H] already discussed in ref. ^1 devoted to 
the direct contribution. This allows us to make a connection between the results 
presented here and those obtained in [15 . A more complete phenomenological study 
of the new HI data P| will be presented in a future paper ^7\, in which we shall 
also discuss in detail the link between the present NLO cross section and the cross 
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Figure 4: The structure function F^^j/a compared to L3 data |lj with = 3.7 GeV^ 
and K"^ = 120 GeV^. Statistical and systematical errors are added linearly. 

section based on the exchange of Reggeized gluons jHHl in the t-channel |16| [T7j 

As in paper [THj, we use the MRST 99 (upper gluon) distributions for the parton 
in the proton jSH] and the KKP fragmentation functions jlU]- The strong coupling 
constant is given by an exact solution of the two loop-renormalization group equation 
and we use = 300 MeV. We take Nf = 4. Our calculations are performed at 
= 300.3 GeV and the forward-vr'' cross section is defined with the following cuts. 
In the laboratory system a 7r° is observed in the forward direction with 5° < O^^o < 
25° ; the laboratory momentum of the pion is constrained by x^o = Et^o/Ep > .01, 
and an extra cut is put on the tt^ transverse momentum in the 7* — p center of mass 
system : E^^o > 2.5 GeV. The inelasticity y = Q^/xb^S is restricted to the range 
.1 < y < .6. We consider only the contribution coming from transversely polarized 
virtual photons, we shall comment briefly on the scalar contribution below. 

Our numerical results obtained for the distribution da/dxg^ measured by HI [Hj 
in the range 4.5 GeV^ < < 15 GeV^ are shown in Fig. 5. In order to shorten 
the numerical calculation we do not integrate over Q^, but instead use the average 



Cf = 3.7 GeV 2 
= 120 GeV 2 




0.2 0.4 0.6 0.8 1 

X 
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value of Q^, < >= 8 GeV^, over the above range. We use the scale + in 
the entire series of calculations and we work in the MS factorization scheme. 
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Figure 5: The cross section da/dxB corresponding to the range 4.5 GeV^ < < 
15 GeV^ compared to HI data |8j. 



The direct HO corrections from which the resolved contribution is subtracted, 
called HOs, are different from those obtained in ref. in which we work in the 
virtual factorization scheme. In both schemes they are very large. In ref. JSj 
we noticed that the largest contribution to these corrections comes from the sub- 
processes 7* + 5' 9 + 1 + Q and 7* + ^ — > q + q + q. The sum of the HO^ 
contributions and of the resolved Born contribution should be factorization scheme 
independent, up to corrections. To check this point, let us consider the bin 

2.9- 10^"^ < Xg^ < 3.9-10"^. In ref. ^3] we used the virtual factorization scheme and 
we obtained da/dzs. = 155.5 nb +52.5 nb = 208.0 nb for the sum. Note that the 
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parton distributions used in that case are simply the lowest order distributions (fT^ 
without QCD evolution. In the MS scheme we have da/dxg. = 125.9 nb +95.9 nb = 
221.8 nb, but we use the NLO parton distributions. One can check that the small dif- 
ference between the two sums comes mainly from the gluon distributions not present 
in the lowest order expression. If we only use the quark MS NLO distributions, we 
obtain da^'^ /dxg. = 83.1 nb for the resolved contribution and da/dxg^ = 209 nb for 
the sum. This result shows that the QCD evolution is negligible in this kinematical 
range (besides the generation of a small gluon distribution) and that expression (jl5|) 
gives a good description, in the virtual factorization scheme, of the parton distribu- 
tions in a virtual photon. A similar observation has been made by the authors of 
ref. ^2 in the case of jet production. Because of this small evolution, we also have 
q^{n, M^, Q^) ^ Qoin). Therefore, it is not necessary to subtract the scalar resolved 
component from the direct term and to introduce a scalar (QCD evolved) resolved 
contribution. 

The next point to observe from Fig. 5 is the importance of the HO resolved correc- 
tions compared to the Born resolved contributions, leading to a ratio NLO / Born 2 
independent of Xg. . These large HO corrections correspond to a small value of -E^4 
due to the small cut-off E^^o > 2.5 GeV. For a larger cut-off, for instance E^^o > 
5 GeV, we obtain NLO /Born = 1.65 in the range 2.9 ■ 10"^ < Xg^ < 3.9 ■ 10"^. 
The total cross section is in good agreement with data, slightly overshooting them 
xb ^ 4 ■ 10~^, and little room appears to be left for a BFKL-type contribution 
|17j . However this last statement depends on the scale used in the calculation, here 
M = = Mp = n = {Q"^ + Ej_^)^ , because the cross section strongly depends 
on the scale fi. In ref. we found that this was due to the importance of the 
subprocesses 7* + 5' g + Q + q and 7* + g q + q + q with a gluon exchanged in the 
t-channel. These processes correspond to the opening of new channels that are not 
present at the Born level. They are of order Ola"^) and sensitive to the value of fi 
since there is no loop contributions at this order to compensate the /^-dependence. 
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However, this remark is not true for the resolved part of these subprocesses as 
soon as HO corrections for the resolved cross section are calculated. For instance 
the subprocess 7* + 5'-^5' + 5' + ^ contains a resolved lowest order contribution 
(7* ^ '¥i) + g^9 + (l + Q] loop corrections to this contribution, corresponding 
to HO corrections to the resolved cross section, generate counter terms in ; 
these in turn compensate the //-dependence of the Born cross section. To check 
this point let us again consider the bin 2.9 ■ 10^^ ^ Xg < 3.9 • 10~^. Keeping 



M = M^ = Mf = {Q'^ + El^)2 fixed and E*^o > 5 GeV, we vary /i = C^Q^ + El^ 
with C ranging from .15 to 1.0. The variations with /i^ of the Born and NLO resolved 
cross section are shown in Fig. 6. 
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Figure 6: The variation with fi = CJQ'^ + of the resolved cross section. 



We see that the behavior of the Born cross section and that of the NLO cross 
section are quite different. The latter has a maximum around C ~ .2 and is more 
stable with respect to the variations of C than the Born contribution. This behavior 
does not occur for the NLO direct contribution which always increases when C 
decreases. Let us also note that this behavior cannot be observed for the cut E] n > 
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2.5 GeV. The HO corrections are too large and we cannot reach a maximum of the 
NLO cross section, even for very small values of C. 

Let us conclude this section by noting another difference with respect to the 
direct term in which the large contributions to the forward cross section come from 
subprocesses involving the exchange of one elementary gluon in the t-channel. In the 
resolved case, the elementary gluon becomes reggeized, due to the HO corrections. 
Therefore, the resolved cross section contains contributions corresponding to the 
exchange of a reggeized gluon in the t-channel. 

6 Conclusion 

In this paper we calculated HO corrections to the resolved part of the DIS cross 
section for the production of large- i^^vr^. This involves the calculations of the HO 
corrections to parton distributions in the virtual photon (of virtuality Q"^) and of 
the HO corrections to the resolved subprocess. 

We discuss the issue of the factorization scheme in detail and we establish the 
inputs of the parton distributions in the MS scheme. Then the NLO parton dis- 
tributions are obtained by solving the DGLAP inhomogeneous evolution equation. 
They are confronted to LEP data. 

Our results for the NLO cross section are compared with HI data for the pro- 
duction of forward large- i?_L7rO . We find a good agreement with data, once the direct 
contribution is added to the resolved one. This result is obtained for renormalization 
and factorization scales equal to + Ej_^o ■ In a study of the scale sensitivity, we 
find that the resolved cross section is less sensitive to the renormalization scale than 
the direct cross section. It is interesting to notice that the authors of ref. [T3j ob- 
tained very similar results in their NLO study of the electroproduction of large- i?^ 
forward jets. 

We conclude that the good agreement between the NLO calculations and data 
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leaves little room, in this kinematical range, for a BFKL-type contribution which 
resums a ladder of reggeized gluon. 
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Appendix 1 

The square of the graphs in Fig. 1 possess single and double poles in k"^. 
Actually, as explained in section 2, the only relevant quantity is the square of graph 
(a) S^'^/k'^. Let us isolate the hard subprocess amplitude /i by using the projector 
defined in ref. ^1 IP = ^^^^ (with r|^' = (1, 0, 0, 1)). On the left-hand side, this 
acts on the hard cross section, and on the right-hand side, on the 'j*qq vertices : 



(eg)' ^ 



,2 , / 2 u2^^V■(l 



+ (m^ - k^ 



4?7 ■ k 

2 (g^r + q-'k^' - 2A;^r) l/ip 



1/^ 



2 



2 1,2 

+ ~ iv^iq" - k") + r/^(g^ - A;^)) \^^\^ (A.l) 
k ■ rj 

with l/ip = Tr{r{^ + ^ + m)r(^ + m)} where F describes the coupling of parton p 
to the quark. All the expressions discussed in section 2 can be obtained from (A.l) 
and the phase space integration. 

I d^k' 5{k'^ -m^) = j d^k 5 ((g - A;)' -rrt')=^j d^ dipkdz 

with the definition z = ^ = and the relation fc^. = —{k'^ — m'^){^ ~ ^) + q^z{l — 
z) — m?. The difference {S^'^ — S'^)/k^ does not lead to singular expressions when 
g2 ^0. 
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